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Q Abstract 

We adapt multilevel, force-directed graph layout techniques to vi- 
sualizing djTiamic graphs in which vertices and edges are added and 
I— I removed in an online fashion (i.e., unpredictably). We maintain mul- 

tiple levels of coarseness using a d5Tiamic, randomized coarsening al- 
gorithm. To ensure the vertices follow smooth trajectories, we employ 
^ djTiamics simulation techniques, treating the vertices as point parti- 

Q cles. We simulate fine and coarse levels of the graph simultaneously, 

' coupling the dynamics of adjacent levels. Projection from coarser to 

I finer levels is adaptive, with the projection determined by an affine 

^ transformation that evolves alongside the graph layouts. The result is a 

Q\ dynamic graph visualizer that quickly and smoothly adapts to changes 

in a graph. 

<N 1 Introduction 

Q Our work is motivated by a need to visualize dynamic graphs, that is, graphs 

from which vertices and edges are being added and removed. Applications 
include visualizing complex algorithms (our initial motivation), ad hoc wire- 
less networks, databases, monitoring distributed systems, realtime perfor- 



X 



d mance profiling, and so forth. Our design concerns are: 

Dl. The system should support online revision of the graph, that is, changes 
to the graph that are not known in advance. Changes made to the 
graph may radically alter its structure. 

D2. The animation should appear smooth. It should be possible to visually 

track vertices as they move, avoiding abrupt changes. 
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D3. Changes made to the graph should appear immediately, and the layout 
should stabilize rapidly after a change. 

D4. The system should produce aesthetically pleasing, good quality lay- 
outs. 

We make two principle contributions: 

1. We adapt multilevel force-directed graph layout algorithms [lWal03ll 
to the problem of dynamic graph layout. 

2. We develop and analyze an efficient algorithm for dynamically main- 
taining the coarser versions of a graph needed for multilevel layout. 

1.1 Force-directed graph layout 

Forced-directed layout uses a physics metaphor to find graph layouts llEad84[ 
|KK89[ |FLM94[ IFR91I1 . Each vertex is treated as a point particle in a space 
(usually or M^). There are many variations on how to translate the 
graph into physics. We make fairly conventional choices, modelling edges 
as springs which pull connected vertices together. Repulsive forces between 
all pairs of vertices act to keep the vertices spread out. 
We use a potential energy V defined b}|^ 

" V ' ' V ' 

spring potential repulsion potential 

where Xi is the position of vertex Vi, is a spring constant, fo is a repulsion 
force constant, and is a small constant used to avoid singularities. 

To minimize the energy of Eqn ([l]), one typically uses 'trust region' meth- 
ods, where the layout is advanced in the general direction of the gradient 
W, but restricting the distance by which vertices may move in each step. 
The maximum move distance is often governed by an adaptive 'tempera- 
ture' parameter as in annealing methods, so that step sizes decrease as the 
iteration converges. 

One challenge in force-directed layout is that the repulsive forces that 
act to evenly space the vertices become weaker as the graph becomes larger. 
This results in large graph layouts converging slowly, a problem addressed 
by multilevel methods. 

^Somewhat confusingly, we use V for potential energy as well as the set of vertices. This 
is for consistency with Lagrangian dynamics. 
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Multilevel graph layout algorithms llWal03[ IKCH02II operate by repeat- 
edly 'coarsening' a large graph to obtain a sequence of graphs Go, Gi, . . . , Gm 
where each G^+i has fewer vertices and edges than Gi, but is structurally 
similar. For a pair (Gj, Gj+i), we refer to Gj as the finer graph and Gj+i as 
the coarser graph. The coarsest graph Gm is laid out using standard force- 
directed layout. This layout is interpolated (projected) to produce an initial 
layout for the finer graph Gm-i- Once the force-directed layout of Gm-i 
converges, it is interpolated to provide an initial layout for Gm-2, and so 
forth. 

1.2 Our approach 

Roughly speaking, we develop a dynamic version of Walshaw's multilevel 
force-directed layout algorithm HWalOSn . 

Because of criterion D3, that changes to the graph appear immediately, 
we focused on approaches in which the optimization process is visualized 
directly, i.e., the vertex positions rendered reflect the current state of the 
energy minimization process. 

A disadvantage of the gradient-following algorithms described above is 
that the layout can repeatedly overshoot a minima of the potential function, 
resulting in zig-zagging. This is unimportant for offline layouts, but can 
result in jerky trajectories if the layout process is being animated. We in- 
stead chose a d3mamics-based approach in which vertices have momentum. 
Damping is used to minimize oscillations. This ensures that vertices follow 
smooth trajectories (criterion D2). 

We use standard dynamics techniques to simultaneously simulate all lev- 
els of coarseness of a graph as one large dynamical system. We couple the 
dynamics of each graph Vi to its coarser version V^+i so that 'advice' about 
layouts can propagate from coarser to finer graphs. 

Our approach entailed two major technical challenges: 

1 . How to maintain coarser versions of the graph as vertices and edges 
are added and removed. 

2. How to couple the dynamics of finer and coarser graphs so that 'layout 
advice' can quickly propagate from coarser to finer graphs. 

We have addressed the first challenge by developing a fully dynamic. Las 
Vegas-style randomized algorithm that requires 0(1) operations per edge 
insertion or removal to maintain a coarser version of a bounded degree 
graph (Section |4|). 
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We address the second challenge by using coarse graph vertices as in- 
ertial reference frames for vertices in the fine graph (Section [Sj . The pro- 
jection from coarser to finer graphs is given dynamics, and evolves simul- 
taneously w^ith the vertex positions, converging to a least-squares fit of the 



coarse graph onto the finer graph (Section 3.2.1 1 . We introduce time dila 



tions between coarser and finer graphs, which reduces the problem of the 



finer graph reacting to cancel motions of the coarser graph (Section 3.2.2 1 . 



1.3 Demonstrations 

Accompanying this paper are the following movies g All movies are real- 
time screen captures on an 8-core Mac. Unless otherwise noted, the movies 
use one core and 4th order Runge-Kutta integration. 

• |evl 049_cube.mov!: Layout of a 10x10x10 cube graph using single- 



level dynamics (Section 3.1 ), 8 cores, and Euler time steps. 



• ||evl049_coarsening. mov; Demonstration of the dynamic coarsening 
algorithm (Section |4]). 

• evl049_twolevel.mov: Two-level dynamics, showing the projection 
dynamics and dynamic coarsening. 

• |evl049_threelevel.mov: Three-level dynamics showing a graph mov- 
ing quickly through assorted configurations. The coarsest graph is 
maintained automatically from modifications the first-level coarsener 
makes to the second graph. 



• 'evl049_compare.mov: Side-by-side comparison of single-level vs. three- 
level dynamics, illustrating the quicker convergence achieved by mul- 
tilevel dynamics. 

• [evl 049_multilevel.mov: Showing quick convergence times using mul- 
tilevel dynamics (4-6 levels) on static graphs being reset to random 
vertex positions. 



evl049_randomgraph.movl Visualization of the emergence of the gi- 



ant component in a random graph (Section [6|). (8 cores, Euler step). 

• |evl049_tree.movt Visualization of rapid insertions into a binary tree 
(8 cores, Euler step). 



^ They may be downloaded fro m |http : //ubiety . uwaterloo . ca/ubigraph/ev/' or by 
following the hyperlinks from Figure [l] 
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(a) evl049_cube 



(b) evl049_coarsening 





(c) evl049_twolevel (d) evl049_threelevel 




(e) evl049_compare (f) evl049_multilevel 




(g) evl049_randomgraph (h) evl049_tree 

Figure 1: Still frames from the demonstration movies accompanying this 
paper. 
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2 Related work 



Graph animation is widely used for exploring and navigating large graphs 
(e.g., [HMMOOJ.) There is a vast amount of work in this area, so we focus 
on systems that aim to animate changing graphs. 

Offline graph animation tools develop an animation from a sequence 
of key frames (layouts of static graphs) . Such systems find layouts for key 
frame graphs, and then interpolate between the key frames in an appropri- 
ate way (e.g. [|FE02in . 

In contrast, our system tackles the problem of online or incremental 
graph layout. We use the term online in the sense of online algorithms: 
one receives a sequence of requests (in our case, changes to a graph), and 
each request must be processed without foreknowledge of future requests. 

The key frame approach can be adapted to address the online problem 
by computing a new key frame each time a request arrives, and then interpo- 
lating to the new key frame. For example, Huang et al [HEW98] developed 
a system for browsing large, partially known graphs, where navigation ac- 
tions add and remove subgraphs. They use force-directed layout for key 
frame graphs, and interpolate between them. 

Another approach is to take an existing graph layout algorithm and in- 
crementalize (or dynamize) it. For example, the Dynagraph system [INWOli 
|EGK*04II uses an incrementalized version of the batch Sugiyama-Tagawa- 
Toda algorithm [ISTT8I1I . 

Our basic approach is to develop an incrementalized version of Wal- 
shaw's multilevel force-directed layout algorithm HWalOSlI . We perform a 
continuously running force-directed layout. When a graph change request 
arrives, the graph is immediately updated. This changes the potential func- 
tion, so that the current layout is no longer an equilibrium position. The 
layout immediately starts to seek a new equilibrium. If the changes come 
rapidly enough, the graph will be in continuous motion. One advantage of 
this approach is that changes to the graph are shown immediately, without 
delaying to compute a key frame. 

3 Layout Dynamics 

We use Lagrangian dynamics to derive the equations of motion for the sim- 
ulation. Lagrangian dynamics is a bit excessive for a simple springs-and- 
repulsion graph layout. However, we have found that convergence times 
are greatly improved by dynamically adapting the interpolation between 
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coarser and finer graphs. For this we use generaHzed forces, which are eas- 
ily managed with a Lagrangian approach. 

As is conventional, we write Xi for the velocity of vertex i, and £Cj for its 
acceleration. We take all masses to be 1, so that velocity and momentum 
are interchangeable. 

In addition to the potential energy V (Eqn Q), we define a kinetic 
energy T. For a single graph, this is simply: 

Roughly speaking, T describes channels through which potential energy 
(layout badness) can be converted to kinetic energy (vertex motion). Ki- 
netic energy is then dissipated through friction, which results in the system 
settling into a local minimum of the potential V. We incorporate friction by 
adding extra terms to the basic equations of motion. 

The equations of motion are obtained from the Euler-Lagrange equation: 



d d d 



L = (3) 



dt dxi dx 

where the quantity L = T — V is the Lagrangian. 
3.1 Single level dynamics 

The coarsest graph has straightforward dynamics. Substituting the defini- 
tions of (l]2) into the Euler-Lagrange equation yields the basic equation of 



motion Xi = Fi for a vertex, where Fi is the net force: 

Xi=Y] -K{Xi - Xj) + V —r^ ^ • -rp- ^ (4) 

J ^ ^ llxj^ — xAly \\Xi — Xj\\ 



spring forces repulsion forces 

We calculate the pairwise repulsion forces in O(nlogn) time (with n = \y\, 
the number of vertices) using the Barnes-Hut algorithm I1BH86I1 . 

The spring and repulsion forces are supplemented by a damping force 
defined by = —dxi where d is a constant. 

Our system optionally adds a 'gravity' force that encourages directed 
edges to point in a specified direction (e.g., down). 



7 



3.2 Two-level dynamics 



We now describe how the dynamics of a graph interacts with its coarser 
version. For notational clarity we write for the position of the coarse ver- 
tex corresponding to vertex i, understanding that each vertex in the coarse 
graph may correspond to multiple vertices in the finer graphF] 




In Walshaw's static multilevel layout algorithm [WalQ3], each vertex Xi 
simply uses as its starting position y^, the position of its coarser version. To 
adapt this idea to a dynamic setting, we begin by defining the position of Xi 
to be j/j plus some displacement Si, i.e.,: 



However, in practice this does not work as well as one might hope, and 
convergence is faster if one performs some scaling from the coarse to fine 
graph, for example 



A challenge in doing this is that the appropriate scaling depends on the char- 
acteristics of the particular graph. Suppose the coarse graph roughly halves 
the number of vertices in the fine graph. If the fine graph is, for example, 
a three-dimensional cube arrangement of vertices with 6-neighbours, then 
the expansion ratio needed will be 2^/^ or about 1.26; a two-dimensional 
square arrangement of vertices needs an expansion ratio of w \/2 or about 
1.41. Since the graph is dynamic, the best expansion ratio can also change 
over time. Moreover, the optimal amount of scaling might be different for 
each axis, and there might be differences in how the fine and coarse graph 
are oriented in space. 

Such considerations led us to consider affine transformations from the 
coarse to fine graph. We use projections of the form 



where a is a linear transformation (a 3x3 matrix) and /3 is a translation. 
The variables {a., (3) are themselves given dynamics, so that the projection 
converges to a least-squares fit of the coarse graph to the fine graph. 

^ The alternative to 'j// would be to write e.g. x'^^y meaning the vertex in level I + 1 
corresponding to vertex i in level I. 




Xi = 6i + yi 



Xi = 5i + ay^ 



Xi = 8i + ay, + /3 



(5) 
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3.2.1 Frame dynamics 



We summarize here the derivation of the time evolution equations for the 
affine transformation (a, /3). Conceptually, we think of the displacements 6i 
as "pulling" on the transformation: if all the displacements are to the right, 
the transformation will evolve to shift the coarse graph to the right; if they 
are all outward, the transformation will expand the projection of the coarse 
graph, and so forth. In this way the finer graph 'pulls' the projection of the 
coarse graph around it as tightly as possible. 

We derive the equations for a and using Lagrangian dynamics. To sim- 
plify the derivation we pretend that both graph layouts are stationary, and 
that the displacements 6i behave like springs between the fine graph and the 
projected coarse graph, acting on a and /3 via 'generalized forces.' By set- 
ting up appropriate potential and kinetic energy terms, the Euler-Lagrange 
equations yield: 

o^ = ^^S,yf + y,Sf (6) 

i 

/3 = ^5:<5. (7) 

i 

To damp oscillations and dissipate energy we introduce damping terms of 
—daOc and —dp[3. 



3.2.2 Time dilation 

We now turn to the equations of motion for vertices in the two-level dy- 
namics. The equations for 5i and 8i are obtained by differentiating Eqn 
©: 

Xi = Si + p + avi (8) 

proj. position 

Xi = Si + + api + aiji (9) 

^ V ' 

proj. velocity 

Xi = di+'$ + ayi + laiji + at/i (10) 

^ V ' 

proj. acceleration 
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Let Fi be the forces acting on the vertex Xi. Substituting Eqn ( [To| ) into iti = 
Fi and rearranging, one obtains an equation of motion for the displacement: 



'0 + ayi + 2a.ifi + cxy^ 

proj. acceleration / 



(11) 



The projected acceleration of the coarse vertex can be interpreted as a 
'pseudoforce' causing the vertex to react against motions of its coarser ver- 



sion. If Eqn pi] ) were used as the equation of motion, the coarse and fine 
graph would evolve independently, with no interaction. (We have used this 
useful property to check correctness of some aspects of our system.) 



The challenge, then, is how to adjust Eqn (11) in some meaningful way 
to couple the finer and coarse graph dynamics. Our solution is based on the 
idea that the coarser graph layout evolves similarly to the finer graph, hut 
on a different time scale: the coarse graph generally converges much more 
quickly. To achieve a good fit between the coarse and fine graph we might 
slow down the evolution of the coarse graph. Conceptually, we try to do 
the opposite, speeding up evolution of the fine graph to achieve a better fit. 
Rewriting Eqn ([8]) to make each variable an explicit function of time, and 
incorporating a time dilation, we obtain 

Xi{t) = Si{t) + Pit) + (12) 

^ V ' 

proj. position 

where (p is a time dilation factor to account for the differing time scales of 
the coarse and fine graph. Carrying this through to the acceleration equation 
5delds the equation of motion 



Si 



( 

'0 + + 2cx(t>y,i + acp'^y^ I (13) 
\^ " ' 

\ proj. acceleration 



If for example the coarser graph layout converged at a rate twice that of the 
finer graph, we might take (j) = \, with the effect that we would discount 
the projected acceleration y^ by a factor of 0^ = ^. In practice we have used 
values of 0.1 < (p < 0.25. Applied across multiple levels of coarse graphs, 
we call this approach multilevel time dilation. 

In addition to the spring and repulsion forces in Fi, we include a drag 



term Ff = —dSi in the forces of Eqn (13 ). 
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3.3 Multilevel dynamics 



To handle multiple levels of coarse graphs, we iterate the two-level dynam- 
ics. The dynamics simulation simultaneously integrates the following equa- 
tions: 

• The equations of motion for the vertices in the coarsest graph, using 



the single-level dynamics of Section 3.1 



The equations for the projection a, (3 between each coarser and finer 



graph pair (Section 3.2.1 1 



• The equations of motion 8i for the displacements of vertices in the 



finer graphs, using the two-level dynamics of Section 3.2 



In our implementation, the equations are integrated using an explicit, 
fourth-order Runge-Kutta method. (We also have a simple Euler-step method, 
which is fast, but not as reliably stable.) 



3.4 Equilibrium positions of the multilevel dynamics 

We prove here that a layout found using the multilevel dynamics is an equi- 
librium position of the potential energy function of Eqn ([l]) . This establishes 
that the multilevel approach does not introduce spurious minima, and can 
be expected to converge to the same layouts as a single-level layout, only 
faster. 



Theorem 3.1. Let {X^X) he an equilibrium position of the two-level dy- 
namics, where X = {61,62, ■■■ ,cx, /3,yi,y2, ■■ ■), and X = X = 0. Then, 
{xi,X2, . . . , Xn) is an equilibrium position of the single-level dynamics, where 
Xi = 6i -\- ay I + P, and the single-level potential gradient W = (Eqn 
vanishes there. 

Proof. Since X = 0, the drag terms vanish from all equations of motion. 



Substituting X = and X = into Eqn ([13]) yields Fi = for each vertex. 
Now consider the single-level dynamics (Section 3.1 ) using the positions Xi 
obtained from Xi = 6i -\- ayi + /3 (Eqn (|5])). From Xi = Fi we have have 
Xi = for each i. The Euler-Lagrange equations for the single level layout 
are (Eqn ([3])): 

d dL d ^ 

-L = 



dt dxi dxi 
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Since = dii = 0, we have --£^L = 0. Using L = T - V and that the 

kinetic energy T does not depend on Xi, we obtain 

for each i. Therefore W = at this point. □ 

This result can be applied inductively over pairs of finer-coarser graphs, 
so that Theorem |3 . 1 1 holds also for multilevel dynamics. 



4 Dynamic coarsening 

As vertices and edges are added to and removed from the base graph G = 
{V, E), our system dynamically maintains the coarser graphs Gi, G2, . . . , Gm- 
Each vertex in a coarse graph may correspond to several vertices in the base 
graph, which is to say, each coarse graph defines a partition of the vertices 
in the base graph. It is useful to describe coarse vertices as subsets of V. 
For convenience we define a finest graph Go isomorphic to G, with vertices 
Vo = {{v} : V £ V}, and edges ^0 = {({^^i}, {^^2}) : (^1,^2) G E}. 
Coarser graphs are obtained by merging vertices via set union. The se- 
quence Vb, Vi, . . . , of of coarse graph vertices forms a chain in the lattice 
of partitions of V; we need to maintain this chain of partitions as changes 
are made to the base graph Vq. 

We have devised an algorithm that efficiently maintains Gj+i in response 
to changes in Gj. By applying this algorithm at each level the entire chain 
Gi, G2, . . . , Gm is maintained. 

We present a fully dynamic. Las Vegas-style randomized graph algorithm 
for maintaining a coarsened version of a graph. For graphs of bounded de- 
gree, this algorithm requires 0(1) operations on average per edge insertion 
or removal. 

Our algorithm is based on the traditional matching approach to coarsen- 
ing developed by Hendrickson and Leland I1HL95I1 . Recall that a matching 
of a graph G = {V, E) is a subset of edges MCE satisfying the restriction 
that no two edges of M share a common vertex. A matching is maximal if 
it is not properly contained in another matching. A maximal matching can 
be found by considering edges one at a time and adding them if they do 
not conflict with an edge already in M. (The problem of finding a maximal 
matching should not be confused with that of finding a maximum cardinality 
matching, a more difficult problem.) 
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4.1 Dynamically maintaining the matching 

We begin by making the matching unique. We do this by fixing a total order 
< on the edges, chosen uniformly at random. (In practice, we compute < 
using a bijective hash function.) To produce a matching we can consider 
each edge in ascending order by <, adding it to the matching if it does not 
conflict with a previously matched edge. If ei < 62, we say that ei has 
priority over 62 for matching. 

Our basic analysis tool is the edge graph G* = {E, S) whose vertices are 
the edges of G, and eiSe2 when the edges share a vertex. A set of edges 
M is a matching on G if and only if Mis an independent set of vertices in 
G*. From G* we can define an edge dependence graph £ = {E, which is 
a directed version of G* : 

ei ^ 62 = (ei < 62) and eiSe2 (share a common vertex) 

Since < is a total order, the edge dependence graph £ is acyclic. Figure [2] 
shows an example. 

Building a matching by considering the edges in order of < is equivalent 
to a simple rule: e is matched if and only if there is no edge e' G M such that 
e' e. We can express this rule as a set of match equations whose solution 
can be maintained by simple change propagation. Let m : E ^ {-L,T} 
indicate whether an edge is matched: m(e) = T if and only if e G M. The 
match equation for an edge e G E is: 

m(e) = -im(e') 

e' ; e'— >e 

where by convention /\ = T. 

To evaluate the match equations we place the edges to be considered for 
matching in a priority queue ordered by <, so that highest priority edges 
are considered first. The match equations can then be evaluated using a 
straightforward change propagation algorithm: 
While the priority queue is not empty: 

1. Retrieve the highest priority edge e = {vi,V2) from the queue and 
evaluate its match equation m{e). 

2. If the solution of the match equation has changed, then: 

(a) If m(e) = T, then match(ui, ^2). 

(b) If m(e) = _L, then unmatch(i;i, i;2). 
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Figure 2: Illustration of dynamic coarsening. Top left: A graph with double 

lines indicating matched edges; Top right: its coarsened version, obtained 
by contracting matched edges; Middle: the edge dependence graph for the 
initial graph, using the (arbitrarily chosen) order eg < 62 < 63 < ey < 
65 < ei < 64. Initially the matching is {e2,ee}. Removing edge 62 (bottom 
left) triggers re-evaluation of the match equations for 65, ei, and 64. The 
new matching is {66,65,61}, resulting in changes to the coarsened graph 
(bottom right) . 
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Both match(e) and unmatch(e) add the dependent edges of e to the queue, 
so that changes ripple through the graph. 

Figure |3] summarizes the basic steps required to maintain the coarser 
graph (y',E') as edges and vertices are added and removed to the finer 
graph. 

4.2 Complexity of the dynamic matching 

The following theorem establishes that for graphs of bounded degree, the 
expected cost of dynamically maintaining the coarsened graph is 0(1) per 
edge inserted or removed in the fine graph. The cost does not depend on 
the number of edges in the graph. 

Theorem 4.1. In a graph of maximum degree d, the randomized complexity 
o/addEdge(e) and removeEdge(e) is O^de^^). 

We first prove a lemma concerning the extent to which updates may 
need to propagate in the edge dependence graph. As usual for randomized 
algorithms, we analyze the complexity using "worst-case average time," i.e., 
the maximum (with respect to choices of edge to add or remove) of the ex- 
pected time (with expectation taken over random priority assignments). For 
reasons that will become clear we define the priority order < by assigning 
to each edge a random real in [0, 1], with 1 being the highest priority. 

Lemma 4.2. Let G* = {E, S) be a graph and p : £' ^ [0, 1] an injective 
function chosen uniformly at random assigning to each vertex a real priority in 
[0, 1]. If G* has maximum degree k, the expected number of vertices reachable 
from any vertex in E following only edges from higher to lower priority vertices 
is O(e^). 

Proof. It is helpful to view the priority assignment p as inducing a linear 
arrangement of the vertices, i.e., we might draw G* by placing its vertices 
on the real line at their priorities. We obtain a directed graph {E, by 
making edges always point toward zero, i.e., from higher to lower priorities 
(cf. Figure |4]) . Note that vertices with low priorities will tend to have high 
indegree and low outdegree. 

We write E[-] for expectation with respect to the random priorities p. 
For e G E, let A'^(e) be the number of vertices of G* reachable from e by 
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Dynamic-Match operations 
We use V to indicate the coarse version of a vertex v. 

• addEdge (e) : Increase the count of {vi,V2) in E' (possibly adding edge 
if not already present). Add e to the queue. 

• deleteEdge(e): If e is in the matching then unmatch(e). Decrease the 
count of {vi, V2) in E'; if this count is then delete this edge from E'. 
If either vi or V2 have no edges except e, then remove vi (resp. V2) 
from the coarse graph. Add all edges e' with e ^ e' to the queue. 

• deleteVertex(v) : for each edge e incident on v, deleteEdge(e). Remove 
V from G. 

• addVertex(v): no action needed. 

• match(e), where e = (^1,^2): For each edge e' where e e' , if e' is 
matched then unmatch(e'). Delete vertices vl and V2 from the coarser 
graph. Create a new vertex vi U V2 in G' . For all edges e = {v,v') in 
G incident on vi or V2 (but not both), add a corresponding edge to or 
from vi U V2 in G' . For each e' such that e e', add e' to the queue. 

• unmatch(e), where e = (^1,^2): delete any edges in G' incident on 
viVJv2- Delete the vertex vi U V2 from G' . Add new vertices vi and V2 
to G' . For each edge incident on vi or V2 in G, add a corresponding 
edge to G'. For each e' such that e ^ e', add e' to the queue. 



Figure 3: Basic operations of the dynamic coarsening algorithm. 




Figure 4: A linear arrangement of the edge dependence graph of Figure |2} 
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following paths that move from higher to lower priority vertices. We bound 
the expected value of A^(e) given its priority p(e) = i]: we can always reach 
e from itself, and we can follow any edges to lower priorities: 

E[A^(e) I p{e) = r/] < 1 + ^ Pr(p(eO < r]) -E[N{e') \ p{e') < rj] (14) 

e' : eSe' ' ' 

Let f{rj) = sup,eijIE[A^(e) I p{e) = r?]. Then, 

E[iV(e') I p{e') <v]< P i]~^f{a)da (15) 
Jo 

where the integration averages / over a uniform distribution on priorities 
[0, 1]). Since the degree of any vertex is < k, there can be at most k terms in 
the summation of Eqn (14). Combining the above, we obtain 

f{rj)=supE[N{e)\p{e)=7]] (16) 
< sup ( 1 + V r]E[N{e') \ p{e') < r?] ) (17) 

V e'Se J 

rn 

<l + krj / r]-^f{a)da (18) 
Jo 

<l + k [ f{a)da (19) 
Jo 

Therefore /(r/) < g{ri), where g is the solution to the integral equation 

g{n) = l + k g{a)da (20) 
Jo 

Isolating the integral and differentiating yields the ODE g{r]) = k~^g'{rj), 
which has the solution g{ri) = e^^, using the boundary condition g{Q) = 1 
obtained from Eqn (20). Since < 77 < 1, g{vi) < e^. Therefore, for every 
e G the number of reachable vertices satisfies E[A^(e)] < e^. □ 

Note that the upper bound of O(e^) vertices reachable depends only on 
the maximum degree, and not on the size of the graph. 
We now prove Theorem |4.1 



Proof. If a graph G = {V,E) has maximum degree d, its edge graph G* 
has maximum degree 2{d — 1). Inserting or removing an edge will caus e us 



to reconsider the matching of < e^^'^ edges on averages by Lemma 4.2 
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If a max heap is used to implement the priority queue, 0(de^") operations 
are needed to insert and remove these edges. Therefore the randomized 
complexity is 0(c?e^'^). □ 

In future work we hope to extend our analysis to show that the entire 
sequence of coarse graphs Gi, G2, . . . , Gm can be efficiently maintained. In 
practice, iterating the algorithm described here appears to work very well. 

5 Implementation 

Our system is implemented in C + +, using OpenGL and pthreads. The graph 
animator runs in a separate thread from the user threads. The basic API is 
simple, with methods newVertex() and newEdge(vi, ^2) handling vertex and 
edge creation, and destructors handling their removal. 

For static graphs, we have so far successfully used up to six levels of 
coarsening, with the coarsened graphs computed in advance. With more 
than six levels we are encountering numerical stability problems that seem 
to be related to the projection dynamics. 

For dynamic graphs we have used three levels (the base graph plus two 
coarser versions), with the third-level graph being maintained from the ac- 
tions of the d5mamic coarsener for the first-level graph. At four levels we 
encounter a subtle bug in our dynamic coarsening implementation we have 
not yet resolved. 

5 . 1 Parallelization 

Our single-level dynamics implementation is parallelized. Each frame en- 
tails two expensive operations: rendering and force calculations. We use the 
Barnes-Hut tree to divide the force calculations evenly among the worker 
threads; this results in good locality of reference, since vertices that interact 
through edge forces or near-field repulsions are often handled by the same 
thread. Rendering is performed in a separate thread, with time step t being 
rendered while step t+5t is being computed. The accompanying animations 
were rendered on an 8-core (2x4) iMac using OpenGL, compiled with g-l- -I- 
at -03. 

Our multilevel dynamics engine is not yet parallelized, so the accompa- 
nying demonstrations of this are rendered on a single core. Parallelizing the 
multilevel dynamics engine remains for future work. 
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6 Applications 



We include with this paper two demonstrations of applications: 

• The emergence of the giant component in a random graph: In Erdos- 
Renyi G{n,p) random graphs on n vertices where each edge is present 
independently with probability p, there are a number of interesting 
phase transitions: when p < the largest connected component is 
almost surely of size 0(logn); when p = it is a.s. of size Q{rp/^), 
and when p > it is a.s. of size Q{n) — the "giant component." 
In this demonstration a large random graph is constructed by preas- 
signing to all (2) edges a probability trigger in [0,1], and then slowly 
raising a probability parameter p{t) from to 1 as the simulation pro- 
gresses, with edges 'turning on' when their trigger is exceeded. 

• Visualization of insertions of random elements into a binary tree, with 
an increasingly rapid rate of insertions. 

In addition, we mention that the graph visualizer was of great use in de- 
bugging itself, particularly in tracking down errors in the dynamic matching 
implementation. 

7 Conclusions 

We have described a novel approach to real-time visualization of dynamic 
graphs. Our approach combines the benefits of multilevel force-directed 
graph layout with the ability to render rapidly changing graphs in real 
time. We have also contributed a novel and efficient method for dynami- 
cally maintaining coarser versions of a graph. 
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